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ABSTRACT 

A cornerstone of AGN unification schemes is the presence of an opticaUy and geomet- 
rically thick dust torus. It provides the obscuration to explain the difference between 
type 1 and type 2 AGN. We investigate the influence of the dust distribution on the 
Eddington limit of the torus. For smooth dust distributions, the Eddingtion limit on 
the dust alone is 5 orders of magnitudes below the limit for electron scattering in a 
fully ionized plasma, while a clumpy dust torus has an Eddington limit slightly larger 
than the classical one. We study the behaviour of a clumpy torus at low and high AGN 
luminosities. For low luminosities of the order of ~ 10'*^ ergs^^, the torus changes its 
characteristics and obscuration becomes insufficient. In the high luminosity regime, 
the clumpy torus can show a behaviour which is consistent with the "receding torus" 
picture. The derived luminosity-dependent fraction of type-2-objects agrees with re- 
cent observational results. Moreover, the luminosity-dependent covering factor in a 
clumpy torus may explain the presence of broad-line AGN with high column densities 
in X-rays. 
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1 INTRODUCTION 

The widely accepted unification scheme for active galac- 
tic nuclei (AGN) proposes that the central accretion disk 
and broad line region (BLR) are surrounded b y a geometri- 
cally thick dusty torus (e.g.. lAntonuccilll993l ). The dust in 
the torus obscures the accretion disk and BLR for lines of 
sight which pass through the torus, while they are visible 
otherwise. The spectral energy distributions (SED) of most 
Quasars and AGN in Seyfert galaxies ha ve a pronounced sec- 
ondary peak in the mid- infrared (e.g., ISanders et al.|[T989l : 
lElvis et al.lll994l ). which is interpreted as thermal emission 
by hot dust in the torus. The dust is heated by the pri- 
mary optical/UV continuum radiation, and the torus ex- 
tends from the dust sublimation radius outwards (|Barvainisl 
[l983)- 

To date, the thermal dust emission from the torus 
has only been spatially resolved by infr ared interferometric 
techniques for th e torus in NG C 106 8 ("Weigelt et alj |2004 
IJaffe et al.ll2004l ) and Circinus ( irristr am ot al. 200^). The 
actual geometry and physical properties (e.g., dynamics and 
dust composition) are, thus, still unknown. In particular, the 
geometrical thickness, which determines the covering (ob- 
scuration) factor, remains a puzzle: Although axisymmetric, 
rotating gas and dust configurations with cooling will form 
thin disks, the torus should keep an aspect ratio H /r ^ 0.5 
for most of the AGN activity phase. 

Besides the geometrical thickness, the dynamics of dust 
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in the torus is strongly affected by radiation pressure from 
the primary AGN luminosity. It has been suggested that 
the (radiation) pressure by starformation inside the torus 
may solve the problem of its geometrical thicknes s (e.g., 
lOhsuga fc Umemural Il999l : IWada fc NormanI |2002| ). In a 
competing scenario, the torus consists of a large number of 
small, self-grav itating, dusty molecular c l ouds which form a 
clumpy torus (jKrolik fc BegelmanI Il988l : iBeckert fc DuschU 
[2004 ). 

In this article, we descibe the consequences of AGN ra- 
diation pressure for the clumpy torus model. In Sect. 2, we 
describe the Eddington limit on gas and dust. In Sect. 3 and 
4, we introduce a lower and an upper luminosity limit for 
the torus, respectively. We summarize our results in Sect. 5. 



2 THE EDDINGTON LIMIT FOR THE TORUS 

In the classical picture, the Eddington limit is defined as 
the state when gravity of the enclosed mass balances the 
radiation pressure from the central source, so that 



Ledd = 47rcGAf (r) • — 



(1) 



Here G is the gravity constant, c is the speed of light, and 
m and a are the mass and the cross section of the particle 
which is exposed to the radiation. For a fully ionized plasma 
arround a black hole, the inverse opacity m/cr — is 
dominated by the proton mass and Thomson scattering of 
electrons. This changes in the region of the torus where dust 
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is mixed with gas. For reference, we use kq = ar/nip — 
0.4cm'^/g for the fully ionized gas. Assuming that gravity 
is dominated by the black hole mass M = Mbh = M7 x 
we obtain the classical Eddington limit 



- (std) 



1.26 X lO'^^ergs"^ • M7 



(2) 



The time-averaged AGN luminosity is expected to scale 
with the mass accretion rate in the accretion disk by L = 
r)MAY)(? ■ In a stationary scenario, the accretion disk itself is 
fueled by mass transported through the torus from galactic 
scales. The mass transport rate through the torus A/xorus 
is related to A/ad via A'/ad — A/xorus — Afoutflow This re- 
lation considers mass loss in an outflow or jet during the 
accretion process from the torus towards the inner accretion 
disk. With t — 1 ~ A/outflow /A/xorus, wc obtaiu 



L = TItMtoiubC 



(3) 



The combination of theory and observation for radiative ef- 
ficient accretion with outflows suggests 0.01 V 0-1 and 
0.1 < r < 1 (e.g.. lEmmering et al.lll993 : 1 Pellet ier fc Pudritj 
1 19921 ). In the following we will use a representative value of 
r?T = 0.05. From Eqn. @, we obtain Afxorus = 0.4 Mq yr"^ x 
L45. Here L45 is the bolometric luminosity in units of 
10''^ ergs-^ 

We will now investigate the accretion properties for a 
dusty medium in the AGN torus considering a smooth dust 
distribution and a clumpy structure, respectively. 



2.1 Smooth dust distribution 

Radiative transfer simulations of AGN dust tori frequently 
use dust which is smoothly dist ri buted (e.g..|Pier fc KrolikI 
1 19921 : iGranato fc Danesd 1 1994 ISchartmann et al.l \20(m . 
The radiation which comes from the central AGN directly 
acts on the dust grains in the torus. The absorption cross sec- 
tion of the dust grains in the optical and UV regime can be 
approximated by their geometrical cross section a = nr^^^^. 
Standard size distributions as sume dust grain size s between 
0.025 ^m and 0.25 ^m (e.g.. iMathis et all [l977l '). Using a 
typical dust grain density of 2 — 3 g/cm^ , we obtain an opac- 



ity KDuat ~ 3 — 40 X 10 HQ. Using kdusi in Eqn. ([T}, the Ed- 
dington luminosity for smoothly distributed dust becomes 



-(smooth) ^ Q p. ^ 1 n40 —1 tit 

^edd -0.3-5x10 ergs -Mr 



(4) 

Due to the large value of k, the Eddingtion luminosity 
decreases by 5 orders of magnitude. As a consequence, 
typical AGN luminosities of lO^^ergs"^ would require a 
black hole mass of 10^^ Mq. This is, however, mconsis- 
tent with observed AfsH /^boi-ratios (e.g.. lKaspi et al.ll2000l : 
I Woo fc Urr-^ 120021 '). As a consequence the dust can not be 
gravitationally bound to the black hole. 

The result with the approximated Koust is consistent 
with recent opacity calculations for a m ore realistic dust 
and gas mixture (e.g., ISemenov et al.ll2003 ). They show that 
dust opacities for UV temperatures, which are dominating 
the AGN accretion disk radiation, are about 4 orders of mag- 
nitude larger than the Thomson opacity kq. In principle, IR 
photons coming from more external regions of the torus act 
as a counter force to the the UV photon pressure from the 
AD. However, the IR opacity of the gas and dust mixture 
is ^ 10^ times smaller than the UV opacity, resulting in 
an insigniflcant effect of IR photons when compared to the 



dominating UV photons. Furthermore, the geometry of the 
torus causes the diffuse IR torus radiation to act more or 
less isotropically on the inner wall of the torus (at least in 
the torus plane), so that a possible IR counter pressure is 
even weakened. Contrary, it is rather expected that the main 
effect of th e IR photon p ressure inside the torus is a vertical 
thickening (|Krolikll2007l ). 

Eqn. (|4]) considers dust grains which are decoupled from 
the gas. However, even for a perfect coupling of gas and 
dust in the torus (no drift of the dust relative to the gas), 
the usually assumed mass ratio of gas to dust of 100 raises 
the limit i^^'d°°"'' only 0.001 of the classical limit in 



Eqn. (O. The limit L 



(smooth) 



is valid for an optically thin 



gas dust mixture, while AGN tori are necessarily optically 
thick, providing self-shilding of most of the torus against 
the AGN radiation. This creates a radiation pressure gradi- 
ent at the inner boundary layer (width ruv ~ 1) of the 
torus. The corresponding outward force on this layer is 



L/L 



(smooth) 



(std) / J- (smooth) 



10 times stronger than 



the gravitational pull of the central black hole. This outward 
force would have to be counterbalanced by an enormous in- 
ward pressure gradient in the inner boundary layer. 



2.2 Clumpy dust torus 

An alter native model to smoothly d istributed dust was pro- 
posed bv lKrolik fc BegelmanI (|l988l ). They argue that most 
of the gas and dust in the torus around an AGN has 
to be arrange d in o pticall y thick, self-gravit a ting c louds . 
IVollmer et all (|2004l ) and iBeckert fc Duschll (|2004l ) pre- 
sented a stationary accretion model for the clumpy torus 
(hereafter: SA model), including relations for torus and dust 
cloud properties. The main idea behind this model is that 
clouds are very compact with a large optical depth in the 
UV due to dust which provides self-shielding against the 
AGN radiation and allows the clouds interior to be cold. 
Dust grains on the directly illuminated sides of the clouds 
are exposed to the AGN and individual dust grains are po- 
tentially accelerated and expeled from the cloud due to the 
radiation pressure. Both magnetic fields and dynami cal fric- 
tion of grains in the gas phase of the cloud (see, e.g., ISpitzeil 
Il978l . Sec. 9) can prevent this and transfer the momentum 
to the gas. In a self-gravitating cloud the radiation pressure 
is therefore received by the whole cloud. Thus, the torus is 
limited by the radiation pressure from the central AGN act- 
ing on the dust clouds instead of single grains, so that we 
can define a cloud opacity Kci — -kR^i/Mci- 

The SA model assumes the clouds to be self-gravitating, 
so that the free-fall time equals the sound crossing time 
Rci/cs- This provides a linear relation between cloud mass 
Mci and iici. These clouds should be stable against tidal 
forces in the gravitational field of the central black hole. 
This requir es R%/r^ ^ 2 Mbh/Mci. When combining both 
limits (see Be ckert fc DuschLi2004i l. one finds an upper limit 
for the cloud size 



Rc 



„3/2 



/8G Afi/2 



and a corresponding mass 



(5) 



(6) 



AGN dust tori at low and high luminosities 3 



Here Cs is the cloud-internal speed of pressure waves which 
is of the order of lkms~^. We use this value as the unit 
for Cs in the following. This speed characterizes the cloud 
internal pressure which is required to balance self-gravity, 
and can be understood as the speed of supersonic turbulence 
in the clouds. Alternatively, the clouds may be magnetically 
supported 0- Due to their large cross section, these clouds 
dominante the absorption, scattering, and IR re-emission. 
From the relations for R^i and M^i, we get an upper envelope 
for the opacity of Kci — O.Tkq ■ rpi^ /(csM^^^). The distance 
from the black hole rpc is measured in pc and the speed 
Ca in km s~^. For clouds smaller than the shear limit, the 
opacity Kci cx Rd becomes smaller. With Eqn. ([l]), we obtain 
the Eddington limit for clouds in a clumpy torus, which are 
directly exposed to the primary AGN radiation. 



-(cl) 



1.78 X 10*Srgs"^ • 



CsM 



3/2 



3/2 



(7) 



This is of the same order as in the classical Eddington limit 
(Eqn. ^ and is consistent with observed AGN luminosities 
and black hole masses. Since I/edd oc k,~^, the Eddington 
limit for small clouds is even larger than in Eqn. ((Tjl. 

Eqn. Q shows that L^^'] oc r"^/^. This implies that at 
larger distances, self-gravitating clouds which are directly 
exposed to the AGN radiation become unbound by the ra- 
diation pressure. Thus, distant clouds have to be shielded 
against the AGN radiation by clouds at small radii. As a 
consequence, there should be no significant vertical flaring 
for a clumpy torus; i.e., we expect H/r « const. Further 
consequences of this behaviour will be discussed in Sect.|4l 



3 THE TORUS AT LOW AGN LUMINOSITIES 

In the previous section, we argued that the Eddington limit 
for clumpy dust tori is well in agreement with the range 
of observed AGN luminosities and black hole masses. To 
be clumpy, a torus requires a small volume filling factor 
$v <C 1 for th e dusty clouds. In the context of a SA model, 
iBeckert fc Duschl ( 2004 ) find a mass transport rate through 
the torus Mxorus = SivvY:, where v = ^^^'2 -ff^^Kcpicr is the 
effective viscosity for a torus and E is the surface density. For 
an obscuring torus, the scale height, H, cannot be smaller 
than the mean free path of clouds H ^ I = {4/3)Rc\/^v- 
Otherwise the torus would become transparent for AGN 
photons. The parameter r in the viscosity prescription mea- 
sures the ratio r = l/H. The geometric thickness of the 
torus and the viscosity is maximised for r = 1. For H ^ / 
the cloud density in the torus growth rapidly and the torus 
would collapse to a thin disk. We therefore adopt H — I 
for a working model. After replac ing the mean free p ath b y 
the appropriate expression from iBeckert fc Duschl (|2004l ) , 
we get $v in terms of Mtofus, 



^ The mass density of clouds from Eqn. l[5j and l[6j is p = 
1.6 X 10~^®g cm~^ M-^rpc ■ A dynamically relevant S-field will 
have a strength of a few mG. This leads to gyration times of 
grains shorter than the dynamical time of the clouds max/cs. 
The large densities n/j ~ 8 X 10^ cm^'^ M^Tp^ are sufficient to 
effectively transfer momentum from grains to the gas and to limit 
the drift velocity of grains to about Cs . Both mechanisms support 
the above claim that radiation pressure acts on the whole cloud. 



<l?v = 



.7/4 



3/2 



1/2 



(8) 



The volume filling factor only depends on the mass transport 
rate through the torus, which we parametrized by Mtoius = 
0.4Mq yr"i xL45-(r?r/0.05)"^ (see Sect.[2l). By substituting 
AfTorus, we obtain a hard lower luminosity for the existence 
of an obscuring torus according to the SA model, 



Liov, = 5 X 10''^ ergs 



rjT 



0.05 



(9) 



at which "I>v ~ 1. For clumpy obscuring tori as described, 
it is necessary that the AGN luminosity is L 2> iiow. If 
L > Liow, the volume filling factor becomes $v 1. At 
this point, the SA model would require that the torus col- 
lapses to a geometrically thin disk. As a consequence, most 
of the dust would be driven away (see Sect. I2.l|l . It is, how- 
ever, known that this situation can be avoide d: Lower lumi- 
nositie s go along with lower accretion rates. IVoUmer et al.l 
(2004') showed that for low mass accretion rate, a clumpy 
and almost transparent (/ ^ H) circumnuclear disk (CND) 
can form similar to what has been found aro und the central 
black hole in our Galaxy (|Gusten et al.|[l987h . The difference 
between the clumpy torus and the CND is that the latter 
one loses most of its obscuration properties while there can 
still be IR reprocessing. 

Several observational studies show that at about 
lO^^ergs"^, the Lboi — ^mir- or Lx — I/MiR-relation show 
a significant ch ange in behaviour compared to h igher lu- 
minosities (e.g.. lLutz et al.ll2004l : iHorst et aLlbOO e'l . Appar- 
ently, the main source of MIR emission at L ~ 10'*^ ergs~^ 
is not the proposed, geometrically thick torus anymore. 

A similar low-luminosity limit has been found for 
models where the dust clouds are not produced in a 
torus but released into a wind from an accretion disk 
( Elitzur fc Shlosman|[200^ '). The cutoff at lower luminosities 
is a result of the fact that the mass outfiow rate in the wind 
cannot exceed the mass acc retion rate in the disk. Ta king 
the same r and r) as used in lElitzur fc Shlosmanl l|2006l) , we 
obtain Liow = 2 x 10*^ ergs^^. 



THE DUST TORUS IN THE HIGH 
LUMINOSITY REGIME 



In Sect. 12.21 we have shown that the Eddington luminos- 

f cl) 

ity for clouds in the clumpy torus, I/g^^, depends on the 
cloud- AGN distance as r~^^'^ . A large fraction of AGN radi- 
ate close to or at the classical Eddi ngton limit for Thomson 
scattering (|McLurefc Dunlodl200ll ). We therefore scale the 
actual luminosity to the classical limit L = i'^'^VEdd, where 
^Edd ^ 1 is the Eddington ratio for the AGN. Once L^^^i 
comes smaller than L, the clouds of corresponding Rdir) 
which are directly exposed to the radiation of the central 
source can no longer resist the radiation pressure. This de- 
fines the condition i/i^dd < ^ ^"^^ existence of dust 
clouds of radius Rci{r) in the AGN radiation field. Since 
^cdd r-dependent, we obtain a maximum distance from 
the AGN, rinax(J?ci), at which the largest clouds can with- 
stand the radiation pressure: 

-2/3 



i(-Rci) = 1.3 pc • 



-2/3^2/3 ^1/3 



7 



R. 



(10) 
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The factor (-Rci/^ci,max) ^ 1 accounts for different cloud 
radii _Rci up to the sheer limit -Rd.max (see Eqn. ([5])). As 
mentioned in Sect. 12.21 clouds at distances r > rmax need to 
be shieled, so that no vertical flaring should occur beyond 
'"max- We want to note that rmax does not refer to an outer 
radius of the torus. It relates the maximum possible size of a 
dust cloud to the radiation pressure from the central AGN. 

A proper scaling of the limiting radius, rmax, is in 
units of the dust sublimation radius which presumably 
sets the inner radius of the torus. The sublimation ra- 
dius is depending on the actual dust chemistry and grain 
sizes. A w ell-referenc e d est imation of the r^ub was intro- 
duced by iBarvaini^ l|l987h . providing rsub = 0.4 pc x 

"1^/2 2 8 1/2 

^45 ^sub isoo'^o.os ' where sublimation of graphite grains of 
radius 0.05 /im and a sublimation temperature of 1500 K is 
assumed. For different grain sizes and chemistry, the subli- 
mation temperature and the temperature exponent change. 
Thus, in the case of silicate grains with similar grain size, 
the sublimation radius is larger by about a factor of about 
3. While the actual chemistry, grain sizes and sublima- 
tion tem peratures around AGN are still a matter of de- 
bate (e.g.. lBarvainislll992l : ISitko et al.lll993l : lKishimoto et all 
bOOTl ). reverb eration measurements o f type 1 AGN support 
rsub oc L^^^ (jSuganuma et al1l2006l ). In the following, we 

1/2 

will use the simplyfied relation r^uh = 0.5 pc x L45 , keeping 
in mind the uncertainty due to the actual dust mixture. 

Finally, this allows for a comparison of rmax with the 
inner boundary of the torus: 

^ = 2.3 . c^^L-J^' ■ f-^) (11) 

'sub Y ^cl,max / 

Interestingly, rmax/rsub is approximately of the order of 
unity. That means that at L > L45, the maximum dust 
cloud size can be limited by the radiation pressure rather 
than the shear limit. Which mechanism dominates for an 
individual AGN depends on its actual L/M. We will, thus, 
distinguish between radiation-limited tori (AGN with higher 
£ or L) and shear-limited tori (AGN with lower £ or L) in 
the following. 

An effect of Eqn. pi|) is a change in obscuration prop- 
erties with higher luminosities. In Sect. |3] we briefly sum- 
marized the results from the SA model which requires H = 
I = 4/3-Rci/$v- Contrary to the lower AGN luminosities, 
Rci is now defined by the radiation limit (Eqn. <\11^ ) instead 
of the shear limit (Eqn. ©). From Rd oc L"^/* (Eqn. (fTTj) ') 

and $v oc M^^J^ « ^'^^'^ (Eqn. (|3}), we obtain H oc L^/*. 
Since the torus is expected to have H/R = const, the thick- 
ness of the torus is determined at the reference distance 
R = rsub oc L^^^ . This results in an average thickness of the 
torus, 

Jf/i? oc L-'/^ (12) 

Thus, we expect to see more type 1 AGN at higher luminosi- 
ties for objects which have radiation-limited clumpy tori. 

This result ca n be interpretecl in the framework of the 
"receding torus" (|Lawrenc ? ll99lD : It has been observed 
that for high-luminosity sources, the observed hydrogen col- 
umn density is lower than in low-luminosity sources (e.g. , 
Ueda et al.ll2 003l; 'Bar ger et allbOOSl : iLa Franca et al.ll2005l : 
Akvlas et al] r2006). This is interpreted as a decrease of the 
covering factor of the torus with luminosity; i.e., H/R ap- 



pears to be anti-correlated with L. Recently, ISimpsonI (|2005l ) 
analysed the luminosity-dependence of the type 1 and type 
2 AGN fraction by combining the results from different sur- 
veys. While the original receding torus (LawrencQil99li ') pre- 
dicts a type-2-fraction /2 oc L'^/^ a better fit was found 
for a situation where the height of the torus depends on lu- 
minosity. From Eqn. (|12|) . we would approximately expect 
/2 oc L~°'^^ for radiation-limited dust tori. This is remark- 
ably close to the correlation /2 oc derived by Simp- 
son. We note, however, that for objects with low ^Edd, the 
L~^''*-dependence should not hold but can even be inverted 
(see Eqn. cloudradius) . This would imply that it depends 
on the actual AGN sample properties if a receding torus is 
observed or not. 

Several authors noted that the UV-to-IR dust extinc- 
tion in AGN is lower than what would be inferred from 
the X-ray column density, in particular referring to AGNs 
showing broad-lines in the optical and significant absorption 
in the X-rays (e.g ., I Wilkes et all |2002| : iPerola et all |2004| : 
iBarger et al.ll2005[ l. In the clumpy torus model, the AGN's 
X-ray emission region is completely obscured by an aver- 
age number. A'', of optically thick clouds along a line of 
sight passing through the torus. Since the size of the X-ray 
source is smaller than the typical cloud radius, Rx < Rd, 
the X-ray column density is nx ~ A'' • rid, where rid de- 
notes the column density of an individual cloud. On the 
other hand, the BLR is only fractionally (and statistically) 
obscured by a number of A'' optically thick clouds, since 
Rblb. > Rci- As a result, the average optical depth in the op- 
tical wavelength range, r y, can be approximated by rv ~ 
( jNatta fc Panagi"al Il984l ). Thus, the inferred optical depth 
from the X-rays, A'^ • rd, overestimates the measured opti- 
cal depth TV ~ A for a more extended emission region like 
the BLR. Here, rd ^ 1 is the optical depth of an indi- 
vidual dust cloud. Furthermore, if N is sufficiently small, 
there exists a high probability that parts of the BLR are 
directly visible to the observer, while the X-ray source is 
still obscured (-Rblr > Rci > Rx)- This effect should be- 
come stronger at higher AGN luminosities when the clumpy 
tori beco me radiation-domin ated, since N oc oc L~^^*. 
Actually. IPerola et al.l (|2004l ') report that ~ 10% of their ob- 
served broad- line AGN show high column densities, all of 
them having X-ray luminosities I/2-iokGV > 10*'* ergs~^. 



5 SUMMARY 

We studied the effect of dust on the Eddington limit in 
the molecular dusty torus of an AGN. While the Edding- 
ton limit for smooth dust distributions is approximately 5 
orders of magnitudes smaller than the classical Eddingtion 
limit for a fully ionized plasma, a clumpy dust torus pro- 
vides a similar Lodd — AfBH-relation as the classical one, and 
is in good agreement with observed luminosities and black 
hole masses. The idea of a clumpy torus is based on self- 
gravitating, optically thick dust clouds which are limited in 
size by the sheer of the gravitational potential of the central 
black hole. In the framework of this model, we were able to 
derive a low-luminosity limit for the existence of an obscur- 
ing clumpy torus, which is of the order of L ~ 10*'^ ergs~*. 
Below this limit, the physical and geometrical properties of 
the torus change significantly. Furthermore, we investigated 
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the behaviour of the clumpy torus at high luminosities. We 
found that the largest clouds in the torus become gravita- 
tionally unbound to the central black hole if the AGN radi- 
ates close to the classical Eddington limit. In such a case, 
the dust clouds in the torus are no longer limited in size by 
the shear of the gravitational potential but by the AGN lu- 
minosity. The effective scale height of the radiation-limited 
tori decreases with luminosity, H/R oc The resulting 

L-dependence for the fraction of type 2 AGN, /2 oc L""'^^ 
is consistent w ith an analysis of several AGN surveys by 
ISimpsonI Hooi). We showed that the clumpy torus can ac- 
count for broad-line AGN with high X-ray column densities, 
and that more such objects should be found at high rather 
than at low luminosities. 
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